Chiral Hadronic Mean Field Model including Quark Degrees of Freedom 
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In an approach inspired by Polyakov loop extended NJL models, we present a nonlinear hadronic 
SU(3) (T-uj mean field model augmented by quark degrees of freedom. By introducing the effective 
Polyakov loop related scalar field $ and an associated effective potential, the model includes all 
known hadronic degrees of freedom at low temperatures and densities as well as a quark phase at 
high temperatures and densities. Hadrons in the model exhibit a finite volume in order to suppress 
baryons at high T and ^ This ensures that the right asymptotic degrees of freedom are attained for 
the description of strongly interacting matter and allows to study the QCD phase diagram in a wide 
range of temperatures and chemical potentials. Therefore, with this model it is possible to study 
the phase transition of chiral restoration and deconfinement. In this paper, the impact of quarks 
on the resulting phase diagram is shown. The results from the chiral model are compared to recent 
data from lattice QCD. 

PACS numbers: 12.38.-t, ll.30.Rd, 14.20.Gk, 25.75.Nq 
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I. INTRODUCTION 

Deducing the properties of strongly interacting QCD 
matter at high temperatures and densities is one of the 
major goals of heavy ion collision experiments, as per- 
formed at the Relativistic Heavy Ion Collider (RHIC), 
the Large Hadron Colhder (LHC), or at the upcoming 
Facility for Antiproton and Ion Research (FAIR) . Exper- 
iments at the LHC and RHIC strongly suggest that at 
sufficiently large beam energies, i.e. high temperatures T 
and small baryochemical potentials /is, a new state of 
matter is created, comparable to a nearly perfect fluid 
with very low viscosity [iHll- Possible phase transitions 
in other regions of the phase diagram with larger chemi- 
cal potentials and baryonic densities arc probed at lower 
beam energies. 

Quantum Chromodynamics (QCD) is the theoreti- 
cal framework for the description of strongly interact- 
ing matter. A great difficulty in dealing with the 
QCD Lagrangian results from the fact that the equa- 
tions of QCD cannot be solved analytically and per- 
turbativc approaches on QCD arc suitable only in the 
high-temperature regime. Nevertheless, there is strong 
evidence for the existence of two phase transitions for 
strongly interacting matter. One is the chiral phase 
transition from a state with spontaneously broken chi- 
ral symmetry at low temperatures to a chirally symmet- 
ric phase with vanishing constituent quark masses [H-Q. 
For this transition, the chiral condensate a = (qq) serves 
as a well-defined order parameter. At vanishing bary- 
ochemical potentials ju = 0, lattice QCD methods can 
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be used to study the chiral transition. In this theoret- 
ical approach, the QCD Lagrangian is treated numeri- 
cally on a space-time grid. Calculations with different 
lattice actions consistently suggest a smooth cross-over 
transition with a "critical" temperature in the region of 
Te « 160 MeV At non-zero chemical potentials, 

standard lattice QCD methods arc not applicable due 
to the fermion sign problem. The fundamental phase 
structure of strongly interacting matter is still unknown. 
However, expansion and reweighting methods exist to ex- 
trapolate results from lattice QCD into the region /i > 
of the phase diagram [l^, [lil, [l9l - [28| . In this region, pre- 
dictions on the order and exact location of the phase 
transition from different lattice groups are not consis- 
tent and depend on the used lattice action and other pa- 
rameters such as the lattice spacing and the quark mass. 
While often the existence of a first-order phase transition 
and a critical endpoint at a certain chemical potential is 
predicted [l^, [H, H^-EI , other recent results suggest a 
smooth cross-over transition in a for the whole T-/i-plane 
under investigation [2^ [s^l • 

The other phase transition in QCD matter is the de- 
confinement transition from hadronic bound quark states 
to a quark-gluon plasma (QGP) state at high tempera- 
tures and baryonic densities [331 with the Polyakov loop 
<& as a suitable order parameter. 

Besides direct numerical approaches solving the QCD 
Lagrangian, a broad range of effective theoretical models 
exist for describing known properties of the phase dia- 
gram and thermod yna mic characteristics of strongly in- 
teracting matter [5l. l33 - l40| . In this paper, such a model is 
presented including all known hadronic degrees of free- 
dom [4l| - |43 | as well as a quark phase for the descrip- 
tion of highly excited matter. With this model, it is 
possible to study the QCD phase diagram with a par- 
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ticular focus on the phase transitions. In our approach, 
the description of the hadronic phase is done by a mean 
field sigma-omega type model a quark phase is imple- 
mented following Polyakov-loop-extended Nambu-Jona- 
Lasinio (PNJL) models [H-il. 

This paper is organized as follows. After this introduc- 
tion, we present our quark hadron model and outline its 
theoretical background in Sec. [Ill Results for the order 
parameters for chiral restoration and the deconfinement 
phase transition, as well as thermodynamic properties of 
our model and the susceptibility coefficients are presented 
in Sec, mil A short summary and concluding remarks are 
given in Sec. IIVI 



II. MODEL APPROACH 

For this study, we use a S'C/(3)-flavor a-uj model with 
a non-linear realization of chiral symmetry; see Refs. [stI . 
l4ll . li^ [56j for a more formal and detailed description 
of the original purely hadronic model. In this effective 
model, the Lagrangian in mean field approximation [ssl , 
[53 takes the form 
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where the first term represents the kinetic energy of the 
hadrons as described in Ref. [l^] in more detail. The in- 
teraction of baryons and quarks with the scalar mesons 
fields fj, ( (attractive interaction, see Eq. ([3])) and the 
vector meson fields w, (f) (repulsive interaction) respec- 
tively, is expressed by 



E 



■ipi [70 (ffit^w" -I- + m*] ipi. (2) 



The summation index i runs over the three lightest quark 
fiavors (u, d, s), the baryonic octet, decuplet, and aU 
heavier baryonic resonance states up to masses of tojv* = 
2.6 GeV, which have at least a three star rating in the 
the Particle Data Group listings [H, . 

The mesonic condensates in our model are considered 
in mean field approximation. They include the scalar 
isoscalar field a ^ (uu+dd) (which is identified as the chi- 
ral order parameter) and its strange equivalent ^ (ss), 
as well as the vector isoscalars lj and 4>. The coupling 
strengths gi^x of quarks and baryons to the scalar fields 
account for their effective mass, which to a large extent 
is generated dynamically 



(3) 



The explicit mass term Srrii is chosen such as to reproduce 
the tabulated vacuum masses, which is the mean mass for 
broad resonance states. The lightest baryons exhibit an 
explicit mass of Srrii = 150 MeV and somewhat larger 
values with increasing vacuum mass of the specific par- 
ticle. The quark masses are set to Srriu.d — 6 MeV and 
Snis = 105 MeV. According to Eq. ([31) the decrease of 



the cr-field at high temperatures and densities leads to 
decreasing baryon masses and thus to the restoration of 
chiral symmetry. The resulting effective masses and the 
connected ground state properties are discussed in detail 
at the beginning of Sec. IIIII 

In the same manner, the effective chemical potential 
for quarks and baryons is given by 



fj-* = fM - 9iu:^ - gi<t>4', 



(4) 



and therefore, its value is generated by the couplings Qi^.cj, 
of quarks and baryons to the vector meson fields. 
The mesonic part of the Lagrangian Eq. ([T]) 
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includes the mass terms and the self interactions of the 
vector mesons 
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as well as the self interaction terms of the scalar mesons 
explicitly given by the terms 
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fc2 ( Y + C) + fc^a^C 
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and an explicit symmetry breaking, defined by the terms 
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The scalar dilaton field x can be identified as the gluon 
condensate. It was introduced in order to ensure scale 
invariance of the Lagrangian and to model QCD trace 
anomaly (42j . In the purely hadronic scenarios without 
considering a quark phase, changes in the x-field are very 
small and thus it is kept fixed at its ground state value 
Xo- This does not apply, in the presence of the quark 
phase. In this case, the quarks couple to the dilaton 
field [cf. Eq. (|20p ] and hence, its value changes with the 
temperature. The chosen parameters of the Lagrangian 
are listed in Tab. [J 

The mesonic potential, which sets the dynamically gen- 
erated meson masses, is defined as Vmeson = — 'Cmeson- 
From this, the masses of the mesons in the model are 
obtained by the second derivative of the potential with 
respect to the particular meson species £,j 



(9) 
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fci 






ki 


3.83x10'' MeV^ 


1.40 


-5.55 


-1.07x10'' MeV 


-0.23 



5 


Xo 


94. 


0.067 


401.93 MeV 


38.5 



Table I. General parameters used in the mesonic part of the 
chiral model Lagrangian Eq. ((5|. The couplings of the parti- 
cles to the meson fields are given in Tab. [Til 



In our model the coupling strengths of the baryons 
and he quarks to the fields are the parameters with the 
biggest impact on the resulting mater properties and the 
phase diagram. The couphngs of the baryonic octet arc 
fixed such as to reproduce the well-known vacuum masses 
and nuclear saturation properties (see Refs. [H, |63| for 
details). For the couplings of the baryonic resonances to 
the fields, a straight forward two parameter ansatz is cho- 
sen [4^. This approach avoids having multiple different 
coupling parameters for every single resonance state or 
SU(3) multiplcts. The resonance couplings are connected 
to those of the nuclcons to the respective fields via the 
relations 



(10) 
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where rs and r„ are the so-called scaling parameters for 
the scalar and vector fields. Ref. [i^ demonstrates the 
impact of baryonic resonance states on the phase dia- 
gram. This was essentially done by varying the vector 
coupling strength (change the value of r„, which effec- 
tively controls the abundance of particles at finite bary- 
ochemical potentials in the model). For reasonable cou- 
plings, a first-order phase transition ceases to exist in 
favor of a smooth cross-over due to the continuous pop- 
ulation of heavier resonances states. 

In this study, the vector coupling parameter is kept 
fixed at = 0.9. Furthermore, in order to ensure a 
smooth cross-over of the order parameter and of the ther- 
modynamic quantities at vanishing chemical potentials, 
the scalar couplings are also fixed at a value close to 
unity w 1. For the non- interacting hadron resonance 
gas (HRG), the couplings of all particles to the fields are 
set to zero and 5mi is kept fixed at the respective vacuum 
mass values [H, . 

The couplings of the quarks to the meson fields are 
controlled by gqa, gqQ, Qqu, and g^^, where q stands for 
the respective quark species (u, d, s). The explicit masses 
of the quarks and the nucleons (5mo and their coupling 
strengths to the fields are listed in Tab. [ill In Sec. IIII Al 
a detailed study of the impact of the quark couplings 
on the phase diagram and the nuclear ground state is 
performed. 

All thermodynamic quantities in our model can be de- 
rived from the grand canonical potential, which is defined 

as 
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5<T 9c 


9uj 


u 
d 

s 


6.0 
6.0 

105.0 


-3.5 0.0 
-3.5 0.0 
0.0 -3.5 


4.0 0.0 
4.0 0.0 
0.0 7.1 


N 


150 


-9.83 1.22 


11.56 0.0 



Table II. Explicit particle masses Sttiq and coupling strengths 
of the quarks and the nucleons gi to the fields. In the non- 
interacting HRG scenario, all couplings are set to zero (cf. 
Fig. [5] for detailed parameter study). 



In the heat bath of hadrons and quarks, the thermal 
contribution to the potential includes those from quarks, 
baryons, and mesons 



(13) 



In explicit form, the thermal energy contributions are 
defined as 



(27r)3 
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i(iJ-(fc)+A'*) 
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for the quarks and antiquarks (i € {u, d, s}), which cou- 
ple to the effective Polyakov loop potential $ and its 
conjugate $ as described below. 
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+ e-T(i^;(fc)-A'*) 
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for baryons and antibaryons with the sum running over 
all baryons i?, and 
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i(ijr(fe)-w) 



(16) 



for the mesons, where in this case the sum runs over all 
mesons. The spin-isospin-dcgeneracy factor is denoted 
with 7ij,; for the respective particle species / and 

j li^) ~ y^fc^ + m*'^j ; for the single particle energies. 

The introduction of the quarks in our model is done as 
in other recent extensions of the original Nambu-Jona- 
Lasinio models 0,|43| which include an effective Polyakov 
loop potential. The Refs. O miH, HI, EM^ give a 
comprehensive review of these PNJL models. 

The scalar Polyakov loop field $ in Eq. (fT4|) is retrieved 
by tracing the constant matrix valued temporal compo- 
nent of the SU(3) color gauge background field 



■Tr 



(17) 



Phenomenologically, this field indicates the phase tran- 
sition from a confined hadronic phase to a deconfincd 
partonic phase and therefore from the hadron resonance 
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To [MeV] 


ao 


ai 


a2 


b 


225 - 235 


3.51 


-2.47 


15.2 


-1.75 



Table III. Parameters of the effective Polyakov loop potential 
Eqs. (|18|) . (|19p and coupling parameters of the model. The 
value of To is adapted according to couplings (see text). 



gas to the quark phase. These dynamics are controlled 
by the effective Polyakov loop potential, adapted from 
Ref. [6^ , in the form 

U = a(T)^<P + b ( lnfl-6#$ , , 

2 ^ ^ \T J L (18) 

+4(1.3 + $3) _3(|.$)2] ^ 

which fulfills the Z(3) center symmetry of pure gauge 
QCD. The temperature dependent parameter 



a(T) = ao + fli 



a2 



(19) 



and the parameters therein are chosen such as to repro- 
duce known features of the deconfinement transition . 
At low temperatures, in the confined phase, $ is zero and 
its value rises with increasing temperatures to $ — > 1 
above the critical Polyakov temperature Tq and with in- 
creasing chemical potentials. The model parameters for 
the Polyakov loop potential arc listed in Tab. IIIII For 
the free parameter of the Polyakov loop Tq, a value of 
235 MeV is chosen when used in the non-interacting sce- 
nario, in which all couplings are set to zero, and 225 MeV, 
when all particles in the model couple to the respective 
meson fields as described above. 

Since the frozen dilaton field X = Xo constrains the 
chiral condensate to finite values for all T and /z, a cou- 
pling of the Polyakov loop potential to the dilaton field 
is introduced. This suppresses a when deconfinement is 
achieved. It is realized via 



X = Xo 



$2 $2 



(20) 



a form similar to the one described in Ref. [6J]. This 
coupling leads to a stronger and more reasonable sup- 
pression of the chiral condensate in the presence of free 
quarks in the last three terms of the scalar meson self 
interaction in Eq. 

Minimizing the grand canonical potential Eq. (|12p for 
a given temperature and baryochemical potential with 
respect to its respective fields, leads to a coupled sys- 
tem of non-linear equations of motion for the fields and 
the particle densities Pij^i- In the following, all thermo- 
dynamic quantities are derived starting with the pres- 
sure p = —d^/dV and the entropy density s = dp/dT 
by using the expression for the internal energy e = 
Ts — pV -f /ii Pi, with the summation index i running 
over all particle species in the model. 



In order to use hadronic degrees of freedom at low tem- 
peratures and densities and quarks at higher tempera- 
tures and densities, we introduce an eigenvolume KLfor 
each particle species z, as it was done in Refs. j55ll65ll66|: 
similar approaches to repulsive volume interactions are 
carried out in Refs. |67H7l| . By doing so, it is ensured 
that with increasing densities at higher T and /i - at 
the latest when quark abundances rise quickly at the de- 
confinement phase transition - an effective suppression 
of hadronic degrees of freedom is established due to ex- 
cluded volume effects, as was shown in Refs. [72| - [75| . An 
eigenvolume for every particle of a species i is introduced 
and expressed by 



1 4 

2 3 



7r(2r)3, 



(21) 



with the parameter r defining the effective radius of the 
hadrons in the model. For the relation between the 
baryon and the meson radii, we assume = tbI^- For 
calculations including the quark phase and assuming full 
couplings to the meson fields, a value = 2.7 fm for the 
baryons and V^j^ for the mesons is chosen (cf. Tab. lIV[) . 
Quarks do not exhibit a finite volume. 

In this way, each single hadron occupies a certain vol- 
ume and thereby reduces the total volume V by V^^- This 
results in a modified chemical potential 



Mi 



P 



(22) 



where denotes the partial pressure of species i. 

The particle, energy, and entropy densities are cor- 
rected by the excluded volume and found to be 



3'===^e,/(l + K>0, 

■i 

s' = ^,s,/(l + K>0- 



(23) 



This definition of p* and the volume corrected thermo- 
dynamic quantities ensures that thermodynamic consis- 
tency holds true in the whole r-/i-plane. That means, 
the relation for the internal energy, as given above, and 
all partial derivatives for determining the thermodynamic 
quantities remain valid. 



III. RESULTS 
A. Model Properties 

First, we present results from the described model at 
vanishing baryochemical potential, /^s = 0. Figure [T] 
shows the order parameter for the chiral phase transition 
normalized to its ground state value cr/tro as a function 
of the temperature. The order parameter for the decon- 
finement transition $ as a function T is shown in Fig. [5) 
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asqtad, N^=8 
asqtad, N^=12 
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X. int. HRG+q 
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Figure 1. (Color online) Normalized order parameter for the 
chiral condensate a jai:! as a function of the temperature com- 
pared to lattice data from Refs. 76-83]. The red solid line 
depicts the results from the model for the interacting hadron 
resonance gas (HRG) without a quark phase. The interacting 
HRG with a quark phase is shown by the blue dashed line 
and the green dashed line shows the results for hadrons that 
do not interact with the mesons fields, together with the de- 
scribed quark phase. Comparing the curves, it clearly shows 
that even when including quarks, the chiral phase transition 
is mainly driven by hadronic degrees of freedom. 



The results from our model are compared to lattice QCD 
results (data points or gray band) using different im- 
proved staggered fermion actions. In detail, these are 
the asqtad action from Ref. [zll, the HISQ action from 
Ref. |Z3,[zi,[il, the p4 action from Ref. [zllli,!!!, and 
the stout action from Ref. j8ll-l83l|. For the calculations 
with the described chiral model, here and in the follow- 
ing, four different scenarios, i.e. basic parameter sets with 
the following abbreviations, are distinguished; 

HRG: describes the pure hadron resonance gas without 
excluded volume suppression and a quark phase 
that docs not interact with the mcsonic fields. 

Int. HRG: same as above with the only difference, that 
all hadrons couple to the meson fields as described 
above. 

HRG+q: includes the hadrons with excluded volume 
effects and a quark phase. All couplings to the 
mesonic field are set to zero only the quarks couple 
to the Polyakov loop. 

Int. HRG+q: . Same as above but here, all couplings 
to the meson fields are set to the finite values given 
in Tab. HI 

A summarizing table of the configurations for these sce- 
narios is given in Tab. IIVI For the full interacting sce- 
nario, a value VJ^x = 2.7 fm"^ is chosen for the baryons. 
This corresponds to a mean baryon radius of ~ 



0.86 fm which is slightly smaller than the measured pro- 
ton charge radius r^, = 0.88 fm (ssj . 

The results for the purely hadronic resonance gas that 
is interacting with the meson fields as described above 
(int. HRG) is depicted by the solid magenta line (color 
online). In this case, at the critical temperature Tc = 
164 McV, which quantitatively agrees well with recent 
lattice data, cr falls off in a small temperature range. 
As already pointed out in Ref. [1^, a slower decrease 
of cr, as it is observed in the lattice data, can be achieved 
by taking into account pionic self interactions j86l-[88}. 
Nevertheless, as predicted by all lattice QCD calcula- 
tions, the transition is a smooth cross-over for both or- 
der parameters. As found in Ref. [i^, the inclusion of all 
known hadronic degrees of freedom in a a-io model re- 
stricts the chiral phase transition to a smooth cross-over 
in the whole T-[i plane. However, the first-order liquid- 
gas phase transition for nuclear matter at the ground 
state up to a critical end-point Tc ~ 15 McV stays present 
in our model. 

The second scenario considered here, is the interacting 
HRG with the described PNJL-like quark phase (blue 
dashed line). The decrease of a is also very strong but 
the slope flattens significantly at the critical temperature 
for the Polyakov loop = 175 MeV as shown in Fig. [2l 
However, in comparison to the pure HRG it shows, that 
the slope of (t(T') does not change up to and there- 
fore that the chiral phase transition is mainly driven by 
hadronic degrees of freedom. 

Since quarks in the model do not have an eigenvolume, 
quark degrees of freedom are strongly preferred to be oc- 
cupied at Tc and above. As of this temperature, the num- 
ber of quarks exceeds the number of baryons and mesons 
in the systems. Their smaller coupling strength to the 
chiral condensate (cf. Tab. |ll| leads to a much slower 
decrease in a. On the other hand, the rapid increase 
of the quark multiplicities at Tc generates a distinctly 
steeper slope of the Polyakov loop parameter $ (b) than 
predicted by lattice QCD data, which is compiled in the 
gray band in Fig. [2l 

However the Polyakov loop $ results reproduce those 
from Ref. [4^ depicted by the red circles. More recent 
PNJL models [1^1 also show a significant overshooting of 
the lattice results in the transition region revealing an 
obvious discrepancy between effective models and first- 
principle calculations on this point. In lattice QCD, there 
is no sharp increase in the deconfinement order parameter 
at Tc. 

As two limiting cases, the order parameters for the 
non-interacting HRG is shown with [green dashed line 
in Fig. [2] and without a quark phase. The cr-parameter 
for the latter two scenarios are shown by the flat line 
at (t/cto = 1. The excluded volume parameters of the 
baryons Vi^x for the four different scenarios considered 
had to be adapted in order to both give a smooth tran- 
sition between the two phases and a reasonable agree- 
ment with lattice data for the thermodynamic quantities. 
These requirements cause a rather large eigenvolume of 
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T [MeV] 

Figure 2. (Color online) The Polyakov loop oder parameter 
"I> for the deconfinement phase transition as a function of T 
at /i = 0. For the two scenarios including the quark phase, 
the non-interacting scenario (green line) and the interacting 
scenario (blue line), $ only shows rather small differences 
below the critical temperature Tc ~ 170 MeV. Our results 
are close to those from Ref. ^1] (red points) from which the 
Polyakov loop potential was adopted. Other parametrizations 
of the Polyakov loop potential (e.g. given in Ref. gil) yield 
a very similar slope of (j>. The gray band covers the lattice 
QCD resuhs for $ given in Refs. j6,-M]- 

the baryons in the fully interacting scenario. The val- 
ues chosen for the different model scenarios considered 
in this study are listed in Tab. IIVI In the two scenarios 
with only a hadronic phase, on the other hand, hadrons 
do not carry an explicit eigenvolume and thus excluded 
volume corrections do not apply. 

Figure [3] shows the total densities of baryons and an- 
tibaryons (green lines), mesons (blue lines), and quarks 
plus antiquarks (black lines) at /i = as a function of 
the temperature, for the two scenarios containing a quark 
phase. The results are shown for the non-interacting ideal 
gas scenario (solid lines) and the fully interacting sce- 
nario with the non-zero couplings to the fields (dashed 
lines). Obviously, the coupling of the particles to the 
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HRG+q 


1.8 
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int. HRG-fq 


2.7 
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Table IV. The four different scenarios considered here include 
the non-interacting HRG and the interacting HRG, both with 
and without a quark phase q. For all scenarios, the values for 
the excluded volume are given in fm^. The couplings of 
the baryons B and the quarks q to the scalar and the vector 
meson fields (ges.v, Pqs.v) are either set to zero when marked 
with a cross, or have the finite value listed in Tab. [ll]if indi- 
cated by a check mark (see text). The first two scenarios are 
purely hadronic and do not include a quark phase. 




T [MeV] 

Figure 3. (Color online) Net densities over the temperature 
cubed pi/T^ of all baryons (green lines), mesons (blue lines) 
and quarks (black lines) as a function of the temperature at 
/is = 0. Solid lines depict results for the ideal HRG and 
dashed lines the interacting HRG both including the quark 
phase. For the eigenvolume of the baryons. Vex = 2.7 fm'^ 
is chosen for the interacting case and Vex = 1.8 fm^ for the 
non-interacting HRG (see text). As expected, since there is 
no repulsive vector interaction for the quarks in the ideal sce- 
nario, the quark density below Tc is much higher than in the 
fully interacting scenario. 



scalar fields, which generates the effective masses as de- 
fined by Eq. ([3|), leads to an increase in the particle mul- 
tiplicities at Tc, in particular of the baryon multiplic- 
ity which, however, drops again shortly due to the large 
baryon volume in favor of a quickly rising quark multi- 
plicity. In the temperature range 0.6 Tc < T < 1.3 Tc, 
there is a phase in which quarks and hadrons coexist, as 
expected for a cross-over transition. At a temperature 
T w 1.3 Tc, baryons vanish from the system. 

The quarks, however, show a stronger dependence on 
their couplings to the repulsive vector fields. If the quarks 
do not couple to the repulsive vector field, showing an un- 
physical behavior they appear even at temperatures far 
below Tc and their total density is much larger than the 
baryon density for all temperatures. Obviously, this be- 
havior has a strong impact on the nuclear ground state 
what will be studied below. If the mesons couple to the 
fields, the effective masses of the pseudoscalar mesons 
scale as mj^gg ~ l/tr [i^. Therefore, at temperatures 
T > Tc their multiplicities decrease significantly faster 
than in the non-interacting case due to an increasing 
mass. 

As defined by Eq. ([3]), the effective masses m* of the 
particles in the chiral model are dynamically generated 
by the coupling to the scalar meson fields. Increasing 
temperature and density decreases a and m* accordingly. 
Therefore, chiral symmetry can be restored. Figure [4] 
shows the effective masses of the baryons and of the 
quarks in the fully interacting scenario as a function of 
the temperature. The lower panel (b) shows the effec- 
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Figure 4. (Color online) Effective masses m* of the quarks 
(u,d — q and s), the baryonic octet (lower panel), and the 
decuplet (upper panel) for the fully interacting scenario in- 
cluding hadrons and quarks as a function of the temperature 
at /xs = 0. All shown particles couple to the fields with 
coupling strengths as listed in Tab. [Ill The change of the ef- 
fective masses largely reflects the slope of the scalar cr-field 
(cf. blue line in Fig. [1} ; therefore the values for m' show a 
rather smooth decrease with increasing temperature. 
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Figure 5. (Color online) Vector densities of up and down 
quarks pq (dashed lines) and of nucleons pjv (solid red line) 
at T = as a function of ps/po for zero (a) and non-zero 
quark vector couplings gq^i (b) and different values of the 
scalar couplings gqax (see parameter list below the plots). 
The gray shaded areas indicate possible coexistence of quarks 
and nucleons. If the quark vector couplings take a finite value 
(b), quarks are predominantly present at ps = 4po and all 
nucleons vanish within the model at ps = 15 po. The solid 
black line in (a) depicts pM for the purely hadronic model 
without quarks. To ensure that no free quarks are present in 
the ground state, both couplings have to be chosen sufficiently 
large. 



tivc masses of the baryon oetet {N, A, S, S), the up and 
down valence quarks, denoted as q, and the strange quark 
s. In the upper panel (a) the effective masses of the low- 
est lying baryon resonances in the decuplet (A, S*, S*, 
51) are shown. Since the decrease of the scalar fields hap- 
pens rather smoothly with increasing temperature, the 
decrease of m* is also relatively slow. Nevertheless the 
steepest decline is still visible at Tc [1, H^. For nucleons, 
the mass at T = 1.5 Tc has dropped down to 0.33 itlq , for 
the constituent up and down quarks to 0.27 TOq''' , and 
to 0.69 TOq for the strange quarks. Most baryons vanish 
at T = 1.3 Tc. Therefore the mean field approximation 
should not be affected by too small baryonic masses and 
possible fluctuations as argued in Refs. fool. |9]|. 

In the hadronic version of this effective model, i.e. 
without the quark phase and excluded volume correc- 
tion terms [U, 113 1 the couplings of the baryon octet to 
the mesonic fields and the mesonic potential were chosen 
in such a way as to reproduce the well-known vacuum 
masses of the baryons and the experimentally determined 



masses of the mesons, as well as the nuclear ground state 
properties (e.g. the correct binding energy per nucleon 
and the nuclear compression modulus), and the asym- 
metry energy. The implementation of the quark phase 
is done in a way which ensures that these requirements 
are still fulfilled. Moreover, one has to consider whether 
quarks already appear in the nuclear ground state. 

First, we consider temperature zero. Figure [5] shows 
the vector densities of up and down quarks 

i—u^d 

together with the density of the nucleons at T = 0, 
where the anti-particle density pj is zero, as a function 
of the baryonic density p b for different sets of the quark 
coupling parameters (I-V, see key in figure) in the fully 
interacting model including hadrons and quarks. All den- 
sities arc given in units of the ground state density po- 
In panel (a) the densities for a vanishing coupling of the 
quarks to the vector meson fields gq^ = and for two dif- 



8 



ferent values for the quark - scalar meson field coupling 
idqtrX = —3.5, blue doted line, and gqa.c = —4.0, ma- 
genta dashed line) are shown. It is obvious that without 
the coupling to the repulsive vector field, in both cases 
the quarks appear at densities smaller than the nuclear 
ground state density at pB = 0.15 fm'^. This specific 
choice of quark coupling parameters may thus only be 
physically correct if the quarks in this region of the phase 
diagram are not considered to be unbound particles. The 
green dashed line in both panels of the figure shows the 
density of the nucleons ppf which differs only slightly for 
the different parameter sets in (a) and (b). As a refer- 
ence, the solid black line in panel (a) depicts the nucleon 
density in the purely hadronic model without a quark 
phase. 

The values for the quark scalar coupling in panel (b) 
are the same as before and the value of the up and down 
quark vector coupling is varied between g^^j = 3.0 and 
gquj = 4.0. It shows, that with a scalar coupling of gqa = 
—4.0 for the up and down quarks as well as gqQ = —4.0 
for the strange quarks and a vector coupling equal or 
larger than gq^ = 3.0 for up and down quarks, all quarks 
arc effectively suppressed below the nuclear ground state 
density as expected. Therefore, with the given choice of 
parameters, which is compatible with the corresponding 
parameters in the baryon sector, the quarks do not affect 
the nuclear ground state. In the case of non-zero quark 
vector couplings (b) , quarks are predominantly abundant 
at = 4po a-nd all nucleons within the model cease to 
exist at ps ~ 15 po. 

With this parameter set, the model reproduces the nu- 
clear saturation density po ~ 0.15 fm'^ and the correct 
binding energy per nucleon E/A{po) « — 16 MeV. At the 
saturation density the cornpression modulus of infinite 
nuclear matter is given by [92l | 



(f E/A 



dkj. 



dpi 



<f E/A 



dp2 



(25) 



P=PO 



Calculations of isoscalar giant monopole and dipole reso- 
nances, i.e. compression modes of nuclei, and related ex- 
perimental data suggest a value of K^o = 240 ± 20 MeV 
for the compression modulus of nuclear matter (see 
Ref. [9^ for a recent comprehensive review on how to 
determine Koa)- The purely hadronic model neglect- 
ing the eigenvolume of the particles, gives a value of 
Kao = 225 MeV. However, if quarks are included and 
the excluded volume effects are taken into account in or- 
der to suppress hadronic degrees of freedom above Tci 
the pressure and E/A as a function of the baryon den- 
sity increases much faster depending on the choice of l^x. 
This leads to a significantly higher compressibility mod- 
ulus, which is Koo ~ 370 MeV for an eigenvolume of the 
baryons of Vex = 1 fm^ and an unreasonably high value 
of = 890 MeV if is set to 2 fm^. 

This rapid increase in the density could be overcome by 
making the eigenvolume of the particles temperature or 
density dependent. Phenomenologically, it seems plausi- 
ble that the eigenvolume of one hadron is not perceived 
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Figure 6. (Color online) Energy density over temperature 
to the fourth power at zero baryochemical potential from 
the presented model. Panel (a) shows results from the 
model without quarks. While the non-interacting HRG (black 
dashed line) rises monotonously, in the interacting scenario 
(solid magenta line) e/T* rises quickly at Tc due to the drop 
in the hadron masses. Panel (b) shows the same quantity for 
the model including a quark phase. In the case without the 
meson field interactions (HRG+q - green line), e/T^ grows 
much faster as in the interacting case (HRG+q - blue line), 
due to the non-existent vector repulsion of quarks. Unlike in a 
Hagedorn resonance gas, all quantities in the model saturate 
at very large temperatures caused by a large but still limited 
number of degrees of freedom. 



by another particle further away at low temperatures and 
densities below the nuclear ground state, where parti- 
cles are statistically distributed relatively far from each 
other in a large volume. Therefore, the eigenvolume ef- 
fect would vanish in this region and it would reach its 
actual full value at higher densities, at which the par- 
ticles in the system arc densely packed and would start 
to overlap. This, on the one hand, would still result in 
the same effective suppression of hadrons at high tem- 
peratures and densities in favor for point-like quarks as 
before. However, the introduction of a density or temper- 
ature dependent eigenvolume needs to be done with care, 
as not to violate the thermodynamic consistency of the 
approach (cf. Eq. [zllij, Isl. Particularly, since 

the transition between the hadronic and the quark phase 
is abrupt in the fully interacting model, this seems rea- 
sonable in order to achieve a much smoother transition. 
Work along this line is in progress. 



B. Thermodynamics 

We will now take a closer look at the thermodynamic 
properties of the model. Figure [H] shows the energy den- 
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sity divided by the temperature to the fourth power, 
e/T'^. Panel (a) depicts the purely hadronic model with- 
out quarks in the non-interacting scenario (black dashed 
line) as well as results from the scenario in which hadrons 
fully interact with the meson fields (solid magenta line). 
The interacting HRG shows a rapid rise of the energy 
at Tc- This is due to the drop in the scalar condensate 
a. Therefore, also the effective hadron masses decrease 
rapidly, which leads to large hadron abundances in the 
absence of any volume based particle repulsion. The in- 
crease of the thermodynamic quantities of the ideal HRG 
is obviously much slower above since the masses are 
not dynamically generated via couplings to the scalar me- 
son fields. The purely hadronic model is considered only 
to be valid at least up to temperatures in the region of 

Panel (b) shows the results from the chiral model in- 
cluding the additional quark phase. Below the critical 
temperature, the energy density in the non-interacting 
scenario (green dashed line) is significantly higher than 
in the non- interacting scenario (blue dashed line). This 
is because in this case, there is no suppression by the 
vector meson fields and quarks are not suppressed at low 
temperatures (cf. Fig. [3|). Furthermore, in this scenario, 
the excluded volume of the hadrons T4x is chosen a bit 
smaller (cf. Tab.|TV|. In the region between T = 140 MeV 
and 170 MeV the interacting scenario shows a dip in 
the energy density. This is caused by the eigenvolume 
suppression of the hadrons and a rather slow increase 
of the quark abundances due to their coupling on the 
repulsive vector meson fields. At higher temperatures 
(T > 250 MeV), results from both model scenarios con- 
verge again, reflecting the fact that from this temper- 
ature on the system in both scenarios is dominated by 
quarks. Unlike in a Hagedorn resonance gas [96l |. given 
the large but still limited number of degrees of freedom 
in our model, thermodynamic quantities do not diverge. 

Figure [7] shows the same quantity from the interact- 
ing HRG with and without quarks in comparison to lat- 
tice QCD results from Refs. [zl [zl-il- Up to T^, the 
results from the purely hadronic model are in quantita- 
tively good agreement with newest lattice data with the 
continuum extrapolated stout action [sTI - fssI . Compared 
to the HRG, the somewhat older lattice data from the 
HotQCD collaboration [7^,[7^[83] clearly underestimates 
the energy density in the hadronic regime (T < Tc)- 
In this figure, the impact of the quark phase and the 
associated excluded volume suppression on the thermo- 
dynamic properties of the model becomes obvious. As 
mentioned above, there is a dip in e/T'^ in the vicin- 
ity of Tc- In this region, the rather large eigenvolume 
of the hadrons leads to a flattening of the baryon den- 
sity as a function of the temperature, while the quark 
density is only just starting to rise sharply and cannot 
compensate the lower hadron contribution (cf. dashed 
lines in Fig. [31). In the transition region, results from 
our model including quarks lie in the region between the 
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Figure 7. (Color online) Same quantity as above, again for 
the interacting HRG with (blue line) and without a quark 
phase (ma gent a line). The results are compared to lattice 
QCD data [Tq, I79l483l ) . Up to Tc, the newer stout continuum 
extrapolated data (black points) is in reasonable agreement 
to the pure interacting HRG. Compared to the pure HRG, the 
scenario including quarks shows a suppression of e/T* in the 
region of Tc. This is caused by a slower increase of the energy 
due to the eigenvolume suppression of hadrons and cannot 
be compensated by the quarks which just start to populate 
the system at these temperatures. Therefore, the model with 
quarks lies in the region between newest stout action results 
and the somewhat older lattice data (asqtad and p4 action) 
which show large deviations from the HRG at T < Tc. 



lattice data from the HotQCD collaboration and newest 
continuum extrapolated Wuppertal-Budapcst results. A 
convergence of more recent and highly anticipated data 
from the HotQCD collaboration, probably with physical 
pion masses and continuum extrapolation, towards our 
full-model results seems possible. 

Figure [8] shows the pressure over the temperature to 
the fourth power as a function of T. Depicted are results 
from the presented model (with line styles as in the pre- 
vious figures) compared to lattice data. Here, the same 
observations as for the energy can be made. The pressure 
in both scenarios without quarks rises monotonously and 
again much faster above Tc in the interacting case due 
to the dropping effective masses. The scenarios includ- 
ing quarks differ again, mainly due to their difference in 
the coupling of particles to the vector fields. The pres- 
sure calculated with the interacting model reproduces the 
slope of the stout action results nicely, only undershoot- 
ing lattice results in a small region around Tc- Again, 
compared to HRG results, older lattice data obviously 
yields too small values for T < 1.2 Tc. 

An important quantity, characterizing the phase tran- 
sition, is the trace of the energy momentum tensor e—3p, 
also known as interaction measure or trace anomaly, 
the basic thermodynamic quantity in lattice QCD from 
which other thermodynamic quantities can be obtained 
via integration. Figure [H] shows this quantity (again di- 
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Figure 8. (Color online) Pressure over temperature to the 
fourth power from the chiral model in comparison to lattice 
QCD results a.t n — 0. The line styles depicting the differ- 
ent model scenarios are as in Fig. [6] As for the energy, the 
lattice data show obvious discrepancies between the different 
lattice actions. The model with presumably correct degrees 
of freedom, the interacting IfRG with the quark phase, shows 
reasonably good agreement to latest lattice data (stout, con- 
tirmum extrapolated). 
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Figure 9. (Color online) Interaction measure, (e — 3p)/r*, as 
a function of the temperature at /is = in comparison to lat- 
tice data obtained with different fermion actions. Line styles 
and colors for the different scenarios as depicted in Fig. [G] 
additional lattice data taken from Ref. [ttI.!?^. Up to T ~ Tc 
the data obtained with the continuum extrapolated stout ac- 
tion agrees well with the HRG with quarks scenario. 



vided by T**) as a function of the temperature calculated 
with the chiral model using the four different scenarios 
as described above. As before, the results of the sce- 
narios, which only include a hadronic phase without a 
volume suppression effect, show a rapid rise at T^. The 
results of the other two scenarios, however, show a max- 



imum slightly above Tc- Again, the ideal non- interacting 
scenario (dashed green line) shows a smoother transi- 
tion between the hadronic phase and the quark phase. 
The fully interacting model (dashed blue line) exhibits a 
slowly increasing interaction measure below Tc, compa- 
rable to the dip in the energy density in this region, fol- 
lowed by a rapid rise due to dropping hadron masses and 
a slower decrease above the critical temperature. Results 
from the latter two scenarios show a slope comparable to 
those from lattice QCD. However, our model does not 
reproduce the results of one lattice action consistently. 
A qualitatively good agreement with lattice results using 
the continuum extrapolated stout action is achieved for 
the ideal hadron and quark model (dashed green line). 

A quantity which connects thermodynamic properties 
of the equation of state with the dynamics of heavy ion 
collisions is the velocity of sound Cg which is deduced 
from thermodynamic quantities via 



dp 



(26) 



The square of the speed of sound of nuclear matter equals 
the change in pressure for a given change in energy along 
an isentropic expansion path. Thus, this important prop- 
erty controls the expansion dynamics of the whole fire- 
ball, also of sound-like fluctuations and density perturba- 
tions within the whole expanding system. The expansion 
dynamics and possible high density shock waves evolving 
from the collision of the two nuclei, in particular in un- 
even sized colliding nuclei, are lik ely t o be reflected in 
final state particle observables [97l - |lOC| . Vice versa, the 
determination of the speed of sound from easily acces- 
sible particle observables may provide in-depth informa- 
tion on the thermodynamic properties of nuclear matter, 
especially in relati on to the phase transition where 
significantly drops 10ll4l07 |. 



The velocity of sound squared as a function of the tem- 
perature c alcul ated with the chiral model and from lat- 
tice QCD |l08t is shown in Fig. [TUl Again, the different 
parameter scenarios described above are considered and, 
as expected, it is observed that the smoother the transi- 
tion of the thermodynamic quantities at Tc, the stiffcr the 
equation of state. This seems plausible since a smooth 
change in effective degrees of freedom at the phase transi- 
tion leads to much smoother course in energy density and 
pressure and thus in a much flatter derivative Eq. (j26p . 

The interacting HRG (solid magenta line) shows the 
smallest at Tc due to the rapid increase of degrees of 
freedom in a small range caused by the drop in m* . When 
looking at the speed of sound for the two fully interact- 
ing scenarios, one also observes a second drop down to 
cj. = 0.1 (or a small dip in cj. if the quark phase is taken 
into account) at temperatures slightly above Tc- This 
local minimum is caused by the strange quark-antiquark 
condensate C which falls off at this point and affects m* 
as described in Eq. For a better understanding, the 
(■-field (0.1C(r)/Co) in the interacting HRG scenario is 
illustrated by the thin dashed gray line. The ideal HRG 
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Figure 10. (Color online) Sound velocity squared at zero bary- 
ochemical potential for the different scenarios as described in 
the caption of Fig. [6] as a function of the temperature. For 
both, the pure HRG and the HRG with quarks, the resulting 
equation of state softens considerably if the interactions of 
the particles with the mesonic fields are taken into account. 
The second dip for the interacting HRG at T « 190 MeV, 
which is also slightly visible in the interacting HRG with the 
quark phase, is caused by a sharp drop of the (;^-field (dashed 
gray line). The lattice data (stout, continuum extrapolated) 
is almost flat up to Tc and and rises with a slope comp arab le 
to the HRG+q scenarios for higher T (taken from Ref. [108| |). 



with the quark phase (green dashed hne) shows only 
a small negative deviation from the ideal HRG (black 
dashed line) at Tc- As mentioned above, this is due to 
smooth changes in p and e at the phase transition (cf. 
Fig.[9|). For higher temperatures the value quickly tends 
toward the ideal gas limit cj. = 1/3, as expected. 

In experiment, a distinctive softening of the EoS, i.e. 
a small in the transition region, strongly slows down 
the flow of matter during the collision process and thus 
stops a deflectio n of s pect ator matter. This, so-called, 
burning log effect |l09l . IllOf , causes an observable reduc- 
tion of directed and elliptic flow of particles produc ed in 
sem i-peripheral collisions, as was argued in Refs. 
Ill5| . This behavior was expe rimentally observed by the 
NA49 collaboration Ref. |ll6l | . At beam energies around 
40A GeV, a collapse of the proton directed and elliptic 
flow was found. We note, that there is no softening of 
the EoS in the lattice data. The sound velocity from 
Ref. |l08j | stays constant up to Tc and rises for higher 
temperatures with a slope comparable to chiral model 
results with included quark phase. 

Given these observations, we conclude that our EoS 
based on the interacting HRG including the quark phase, 
offers not only feasible degrees of freedom in a wide range 
of temperatures and densities but also exhibits properties 
consistent with experimental observations. 



IV. SUMMARY 

In this paper, results from a well tested chiral ef- 
fective a-uj model were presented, including all known 
hadronic degrees of freedom, which for this study was ex- 
tended by a quark phase. The quarks were added in an 
PNJL-likc approach in which they couple to an effective 
Polyakov loop potential. The chiral order parameter and 
the Polyakov loop order parameter for the deconfinement 
transition show a critical temperature T « 165 MeV, 
which is in line with all latest lattice QCD predictions. It 
shows, that the chiral phase transition is mainly driven in 
the hadronic regime and there are no qualitative changes 
concerning the chiral transition when adding a quark 
phase. 

At higher temperatures and densities hadronic de- 
grees of freedom are effectively suppressed in favor of 
a purely quark-dominated QGP phase by cigcnvolume 
corrections. Therefore, the presented model offers a de- 
scription of nuclear matter with the correct degrees of 
freedom in a wide range of temperature and baryochcm- 
ical potential. Additionally, with the chiral model it 
is possible to study both, the chiral and the deconfine- 
ment phase transition. Fixing the value for the coupling 
strength of quarks to the repulsive meson vector fields, 
which has a major influence on the structure and the 
properties of the resulting phase diagram, reasonable nu- 
clear ground state properties and a first-order liquid-gas 
phase transition can be reproduced. Nevertheless, the 
value of the nuclear matter compression modulus may 
restrict our choice of values for the finite particle vol- 
ume and the repulsive quark vector couplings. For a 
future study, we intend to introduce a thermodynami- 
cally consistent implementation of temperature or den- 
sity dependent eigenvolumc corrections that could rem- 
edy this problem. For both, the chiral and the deconfine- 
ment phase transition, a smooth cross-over transition at 
vanishing and also finite baryochemical potentials is ob- 
tained, as previously shown without quarks for the chiral 
transition in Ref. [43| . 

The thermodynamic quantities from the chiral model, 
generally show a qualitatively reasonable agreement with 
lattice data, even though, in particular for the interac- 
tion measure there are still significant differences between 
the various lattice QCD results depending on the lattice 
action used. The softening of the sound velocity in the 
transition region was present in all scenarios except the 
pure HRG, where no rapid change of degrees of freedom 
takes places. For the HRG with the included quark phase 
the largest drop in was observed. If all particles in the 
model fully couple to the mesonic fields, the sound ve- 
locity squared gets close to zero at the phase transition. 
This significant softening meets a major constraint on 
the EoS as it is suggested by experimental observations 
on the drop in the fiow by the NA49 collaboration jll6l | . 

The chiral effective model presented here, is able to 
provide an EoS to be used in hydrodynamic models for 
studying the dynamics of heavy ion collisions. This is a 
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major benefit of an integrated model which has the right 
degrees of freedom in a wide range of the phase diagram. 
Therefore, the resulting EoS is well s uited for calcula- 
tions with hybrid models [qqI. [TT7U12H . like the UrQMD 
hybrid model |l22| . since it exhibits known features of 
QCD and has exactly the same degrees of freedom in 
the hydrodynamic stage as in the early and subsequent 
cascade stage. 
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